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Abstract.  By  integrating  the  truncated  complex  scalar  gravitational 
motion  equations  for  an  anelastic,  rotating,  slightly  elliptical  Earth,  the 
complex  frequency  dependent  Earth  transfer  functions  are  computed  di¬ 
rectly.  Unlike  the  conventional  method,  the  effects  of  both  oceanic  loads 
and  tidal  currents  are  included  via  outer  surface  boundary  conditions,  all 
of  which  are  expanded  to  second  order  in  ellipticity.  A  modified  ellipticity 
profile  in  second  order  accuracy  for  the  non-hydrostatic  Earth  is  obtained 
from  Clairaut’s  equation  and  the  PREM  Earth  model  by  adjusting  both 
the  ellipticity  of  the  core-mantle  boundary  and  the  global  dynamical  el¬ 
lipticity  to  modern  observations.  The  effects  of  different  Earth  models, 
anelastic  models,  and  ocean  models  are  computed  and  compared.  The 
atmospheric  contributions  to  prograde  annual,  retrograde  annual  and  ret¬ 
rograde  semiannual  nutation  are  also  included  as  oceanic  effects.  Finally, 
a  complete  new  nutation  series  of  more  than  340  periods,  including  in- 
phase  and  out-of-phase  parts  of  longitude  and  obliquity  terms,  for  a  more 
realistic  Earth,  is  obtained  and  compared  with  other  available  nutation 
series  and  observations. 


1.  Introduction 

Usually,  the  methodology  of  the  nonrigid-Earth  nutation  study  fall  into  four 
categories:  (1)  completely  fitted  from  modern  observation  (VLBI  plus  LLR) 
(McCarthy,  1996);  (2)  semi- analytical  method  (Mathews  et  al,  1991,  2000); 
(3)  based  on  the  Hamiltonian  theory  and  variational  principles,  extending  the 
method  of  rigid  Earth  nutation  study  to  the  study  of  nonrigid  Earth  nutation 
(Getino  &  Ferrandiz,  1999);  and  (4)  purely  geophysical  approaches,  in  which 
the  so-called  nonrigid  Earth  nutation  transfer  function  is  integrated  from  the 
equations  of  Earth  motion  by  a  normal  mode  expansion  method  (Wahr  1981, 
Dehant  &  Defraigne  1997)  or  by  a  direct  numerical  integration  method  (Schastok 
1997,  Huang  1999),  which  is  used  in  this  study  again. 


2.  The  numerical  integration  method 

Smith  (1974)  developed  an  algorithm,  in  which  the  scalar,  but  coupled  due  to  the 
ellipticity,  ordinary  equations,  as  well  as  a  set  of  boundary  conditions,  governing 
the  elastic-gravitational  normal  modes  of  a  rotating,  slightly  elliptical  Earth  in 
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hydro-static  equilibrium,  having  an  isotropic  elastic  constitutive  relation,  are 
finally  derived,  truncated  to  Unite  order  and  described  in  functional  forms,  on 
which  are  based  almost  all  the  modern  nutation  and  tides  theories  for  a  nonrigid 
Earth. 

The  given  equations  of  motion  are  in  ordinary  differential  form.  In  order 
to  overcome  the  resonance  “pollution”  and  to  detect  small  signals,  such  as  the 
contributions  of  atmosphere  and  oceanic  current,  Schastok  (1997)  and  Huang 
(1999)  re-derived  the  boundary  conditions  about  the  stress  tensor  and  gravi¬ 
tational  potential  held  to  second  order,  instead  of  first  order  given  by  Smith 
(1974),  in  ellipticity,  as  well  as  the  boundary  conditions  about  the  displacement 
held  (Huang,  1999;  Huang,  2000a). 

2.1.  The  contributions  of  the  ocean  and  atmosphere 

Unlike  the  conventional  method,  the  effects  of  atmospheric  relative  motions, 
as  well  as  the  oceanic  load  and  current  effects,  are  included  via  non-free  outer 
surface  boundary  condition  in  this  paper. 

If  the  Earth  does  not  include  ocean  and  atmosphere,  the  outer  surface  of 
the  solid  Earth  is  a  free  one,  i.e.  the  pressure  and  the  stress  tensor  at  the 
outer  boundary  of  the  solid  Earth  cancel  in  the  integration.  On  the  contrary, 
if  the  ocean  and/or  atmosphere  are  included  in  the  Earth  system,  the  outer 
surface  of  the  solid  Earth  is  no  longer  a  free  one.  The  stress  tensor  at  the 
outer  boundary  of  the  solid  Earth  is  given  by  Sasao  &  Wahr  (1981)  in  tensor 
format  and  expanded  into  scalar  format  of  three  components  (P™,  QJ1,  RJ1) 
to  second  order  in  ellipticity  by  Schastok  (1997)  and  Huang  (1999).  In  their 
expressions  of  the  continuity  quantity  at  the  non-free  outer  surface  boundary, 
there  are  two  quantities  related  to  the  atmospheric  (or  oceanic)  contribution, 
hjl(a,u)  and  AJl(a,u),  which  represent  the  relative  angular  momentum  and  the 
load  respectively,  where  a  and  u  denote  the  mean  radius  of  the  Earth  and  the 
circular  frequency  of  the  nutation,  and  /  and  m  are  the  degree  and  order  of 
GSSH. 

We  adopt  the  model  of  Zahel  (1995),  which  gives  both  the  amplitudes  and 
the  phases  of  several  dominant  diurnal  tides  of  the  load  and  the  relative  mo¬ 
tion  separately,  and  an  empirical  function,  linear  in  frequency,  to  interpolate 
and  extrapolate  the  other  diurnal  tides  as  Schastok  (1997)  did.  However,  the 
atmospheric  spectrogram  is  usually  contaminated  strongly  and  there  is  not  any 
atmospheric  model  which  can  give  the  amplitude  and  the  phase  of  the  corre¬ 
sponding  diurnal  atmospheric  tides. 

In  this  paper,  we  adopt  the  concept  of  celestial  effective  angular  momentum 
(CEAM)  to  discuss  in  the  frequency  domain.  The  EAM  in  TRF,  y,  can  be 
converted  pointwise  to  the  CRF,  y',  by  multiplying  by  a  factor  of  the  sidereal 
rotation,  (or  etQt)  (Brzezihski,  1994). 

The  original  near-diurnal  signals  in  the  terrestrial  reference  frame  (TRF) 
become  low-frequency  signals  in  the  celestial  reference  frame  (CRF);  on  the 
contrary,  the  original  long-period  signals  in  the  TRF  become  high-frequency 
signals  in  the  space.  Then  we  can  use  a  low-pass  filter  to  keep  the  low-frequency 
signals  of  y'  and  filter  out  the  high-frequency  part  easily  (Bizouard  et  al. ,  1998). 
Huang  et  al.  (2000b)  processed  the  atmospheric  EAM  data  from  1980.0  to 
2000.0  provided  by  NCEP  in  which  there  are  four  evenly  distributed  points 
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Figure  1.  Ellipticity  profile  of  the  Earth  for  equipotential  surfaces. 


everyday.  That  study  assumes  that  the  oceanic  reaction  to  the  atmosphere  is 
non-inverted-barometer  (NIB)  because  it  is  more  true  for  the  high  frequency 
variation  likes  diurnal  tides.  The  amplitudes  of  the  long-trend  term  and  three 
periods  from  standard  least-squares  fit  as  well  as  the  corresponding  relative 
angular  momentum,  /?.,  caused  by  the  wind  are  given  in  that  paper  and  are  used 
in  this  study. 

2.2.  Inelastic  mantle 

If  the  mantle  inelasticity  is  considered,  the  Lame  coefficients  become  frequency 
dependent  and  complex,  and  will  produce  out-of-pliase  contributions  to  nutation 
(Dehant,  1987).  A  frequency-to-the-power-Q'  model  (Anderson  &  Minster,  1979) 
is  used  to  compute  them  in  this  study.  The  radially  dependent  quality  factor, 
Q,  is  taken  from  the  anelasticity  model  QMU  of  Sailor  &  Dziewonski  (1978), 
and  the  reference  frequency  ujo  and  the  power  a  are  taken  as  200s  and  0.15 
respectively  (Wafir  &  Bergen,  1986). 

2.3.  The  ellipticity  profiles  of  the  Earth  interior 

In  addition  to  the  density  profile,  the  ellipticity  profiles  of  the  Earth’s  interior  are 
very  important  for  the  nutation  computation,  although  they  are  less  well-known. 

The  ellipticity  profiles  of  the  Earth’s  interior  depend  on  the  density  profile 
and  can  be  obtained  from  the  second  order  Clairaut’s  equation  for  the  Earth  in 
hydrostatic  equilibrium,  which  is  plotted  as  a  solid  line  in  Figure  1. 

In  order  to  be  consistent  with  modern  observations,  we  assume  that  the 
CMB  and  two  discontinuities  (r  =  5701  and  5971  km)  are  not  in  hydrostatic 
equilibrium  using  two  constraints:  (1)  the  resultant  H  is  consistent  with  the 
observed  value,  i.e.  Hobs  ss  0.00327379  while  from  PREM  (Dziewonski  &  An¬ 
derson,  1981),  JJPREM  ~  0.003240,  and  (2)  the  resultant  retrograde  FCN  period 
(real  part)  is  -432.8  sidereal  days  (with  ocean),  while  it  is  about  -458  days  for 
PREM.  In  order  to  satisfy  this  requirement,  €qmb  is  modified  from  0.002547 
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(for  PREM)  to  0.002666,  while  the  65701  is  modified  to  0.003200  and  65971  is 
modified  to  0.003224.  The  modified  ellipticity  profile  of  the  equipotential  sur¬ 
faces  can  then  be  obtained  from  Clairaut’s  equation  and  plotted  in  Figure  1  also 
(dotted  line). 

As  a  result,  the  period  of  the  prograde  free  inner  core  nutation  (FICN)  is 
also  obtained,  which  is  +465.6  sidereal  days  (with  ocean)  for  modified  ellipticity 
profile  and  is  +470.6  sidereal  days  for  PREM  model. 


3.  Results  and  comparison 

By  integrating  the  truncated  complex  scalar  gravitational  motion  equations,  the 
complex  Earth  transfer  functions  for  each  nutation  frequency  are  computed  di¬ 
rectly  and  are  then  convolved  with  the  rigid-Earth  nutation  series  REN2000 
(Souchay  et  al. ,  1999).  In  the  integration,  the  density  profile  of  the  Earth 
model  PREM  is  used,  and  the  normal  step  length  of  the  integration  is  5km. 
Finally,  the  resulting  nonrigid-Earth  nutation  series  of  343  terms,  in  which 
the  effects  of  the  ocean  load  and  current  and  atmospheric  winds  are  included, 
are  calculated  and  available  for  interested  readers  by  anonymous  FTP  server 
(ftp : / / center . shao . ac . cn/clhuang/nut at ion/hj 12000-343 .txt)  or  by  con¬ 
tacting  the  authors,  because  it  is  too  long  to  present  in  this  paper. 

Moreover,  the  main  four  nutation  terms  (18.66  years,  1  year,  1/2  year  and 
13.66  days)  have  been  compared  with  other  available  observational  and  theoret¬ 
ical  series  (see  Table  1).  They  are:  (1)  IERS  Conventions  (1996)  (McCarthy, 
1996)  nutation  series  (abbreviated  as  IERS96  hereafter)  fitted  from  VLBI  ob¬ 
servations  by  Dr.  T.  Herring  ;  (2)  Dehant  &  Defraigne  (1997)  model  (DD97); 
(3)  Schastok  (1997)  model  (S97);  (4)  Getino  &  Ferrandiz  (1999)  model  (GF99); 
and  (5)  Mathews  et  al.  (2000)  model  (MBH2000). 
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Table  1.  Comparison  of  the  four  main  nutation  terms  (units:  mas) 


fi'tpop 

$£op 

IERS96 

18.66yr 

-17206.277  3.645 

1.553 

9205.356 

S97 

-17206.35 

2.28 

1.11 

9205.05 

DD97 

-17205.39 

— 

— 

9204.87 

GF99 

-17206.251 

3.235 

1.450 

9205.403 

MBH2000 

-17206.458 

3.334 

1.533 

9205.175 

This  paper 

-17206.271 

3.267 

1.309 

9205.147 

IERS96 

147.538 

lyr 

1.121 

-0.198 

7.388 

S97 

147.66 

-0.24 

+0.14 

7.27 

DD97 

147.52 

— 

— 

7.40 

GF99 

147.550 

1.247 

-0.445 

7.256 

MBH2000 

147.557 

1.197 

-0.203 

7.400 

This  paper 

147.567 

0.250 

-0.173 

7.423 

IERS96 

0.5  yr 

-1317.014  -1.400 

-0.464 

573.058 

S97 

-1316.95 

-1.42 

-0.43 

572.95 

DD97 

-1316.94 

— 

— 

572.95 

GF99 

-1317.028 

-1.749 

-0.625 

573.005 

MBH2000 

-1317.159 

-1.382 

-0.459 

573.061 

This  paper 

-1316.971 

-1.435 

-0.432 

572.904 

IERS96 

13.66d 

-227.720  0.269 

0.136 

97.864 

S97 

-227.71 

0.37 

0.18 

97.86 

DD97 

-227.82 

— 

— 

97.90 

GF99 

-227.733 

-0.407 

-0.153 

97.869 

MBH2000 

-227.588 

0.333 

0.163 

97.814 

This  paper 

-227.718 

0.328 

0.178 

97.858 
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